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ABSTRACT 



Fiber light guides are used in conjunction with holo- 
graphic interferometry. A Michelson interferometer was 
constructed to test for temporal coherence of light passing 
through fiber optics. A finite-fringe hologram of a no-flow 
condition was taken to prove the feasibility of applying 
fiber optics to this field of flow investigation. 

The effect of opaque bodies on interferometric data 
inversion was sought. Two schemes for replacing the missing 
fringe data were used to see if the accuracy of the calcu- 
lated density fields could be enhanced. 
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I. EXPERIMENTAL WORK ON FIBER OPTICS 



A. INTRODUCTION 

Holographic interferometry is an invaluable and frequently 
used tool in modern flow field visualization. Numerous 
articles have reported on the application of holography to 
interferometry, among the first of which was presented by 
Heflinger, et al. [Ref. 1], The history, the theory, and 
diverse applications of holography are presented in Collier 
[Ref. 2]. Holographic application to interferometry has been 
extensively investigated for pure gas flows [Ref. 3]> for 
flows around transparent bodies [Ref. 4], and for flows 
around opaque bodies [Ref. 5]* 

The principles of holography have been known since 19^7 
when Gabor developed the principle as an improvement for 
electron microscopy [Ref. 2], but the current scientific 
applications awaited the development of a coherent visual 
radiation source. Previously, flow field visualization 
relied on interferometric techniques utilizing light sources 
of poor temporal coherence. 

The specialized area of interferometric holography 
includes such techniques as real-time, time-averaging, 
stroboscopic, and time-lapse. They are discussed in Ref. 6. 
The time-lapse technique, also referred to as double 
exposure, is used for flow field visualization. This tech- 
nique involves superimposing the phase and amplitude 
information of a flo w condition upon that of a no-flow 



condition. An improvement to time-lapse holography was 
achieved by the use of diffusing screens positioned in the 
scene beam. The idea was advanced by Heflinger [Ref. 1], 
while Gates further developed the concept [Ref. 7]. In 
normal side-band holography, each point on the object acts 
as a diffusely reflecting light source. Light reflected 
from all points on the object can mutually interfere with 
light from all other points. The total interference 
pattern produces a three-dimensional image. This three- 
dimensionality can be obtained in holographic interferograms 
of transparent objects by introducing a diffusing glass 
plate in the scene beam just prior to the test section. 

Two varieties of time-lapse holography are the infinite 
fringe and the finite fringe methods. Infinite fringe 
patterns occur when the wave fronts of the separate expo- 
sures see exactly the same optical arrangement. Finite 
fringe patterns occur when a finite displacement is intro- 
duced In the optical arrangement between exposures. The 
fringe for the undisturbed region are parallel, equispaced 
lines which will be perpendicular to the direction of 
translation. These parallel, equispaced lines form a 
convenient scale for determining the fringe shift values 
in the disturbed region of flow. 

The double exposure has the advantage of relaxing the 
tolerances of. the optical equipment. Both flow and no-flow 
beams travel common optical paths which effectively elimi- 
nates the interference fringes from optical imperfections 
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of the experimental apparatus. This phenomenon is such that 
good quality holograms can be constructed with extremely 
rudimentary setups and optical components. 

One final improvement in technique was required for the 
true versatility of holographic interferometry to emerge. 
Mechanical stability of the experimental components was 
required to achieve proper fringe interference patterns. 

The stability needed to be on an order of the wavelength 
of the light being used. The pulsed laser has eliminated 
stability problems during exposures by its ability to 
"freeze" the flow through exposure times in the tens of 
nanoseconds. Inherent in the pulsed laser is a loss of 
coherency length as compared to the continuous-wave laser 
(c.w. laser). Coherency lengths are strongly dependent on 
the number of modes that can oscillate simultaneously. For 
pulsed lasers, coherency lengths are normally on the order 
of tens of centimeters. 

There is a need to relax the typical experimental 
arrangement to further expand the versatility of holographic 
interferometry. Data for flow field visualization must 
normally be collected from wind tunnel tests which only 
simulate real life. Freedom from the laboratory environment 
could be achieved by adapting fiber light guides to channel 
coherent light to and from real-life flows. Interferometric 
data could thus be collected on many interior flows or 
hazardous flows which do not lend themselves to standard 
experimental setups. Real-life flows such as inlet flows 
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or compressor flows could be investigated. To determine 
the feasibility of applying fiber optics in this fashion, 
a basic interferometric experiment was conducted using a 
fiber light guide to transmit a portion of the laser beam 
to the photographic emulsion. 

B. EXPERIMENT 

1. Equipment 

The initial endeavor to achieve holograms using 
fiber optics was made with a 15 m.w. continuous wave 
helium-neon laser. Having found this inadequate, a Korad 
K-l pulsed ruby laser with a Pockels cell Q-switching 
device was used in its place. The laser operates at a 
wavelength of 6943 Angstroms and has an effective exposure 
time of 20 nanoseconds. The fiber optics employed were 
1/4-inch bundles, one being twelve inches in length and the 
other 24 inches. The c.w. laser was used for alignment of 
the pulsed laser beam through all the optical components. 
Also, the c.w. laser was used to establish the relative 
intensities of the scene and reference beams. Optical 
components consisted of several concave lenses, a beam 
splitter which was 20$ reflective, front surface mirrors, 
and diffusing glasses of several different densities. A 
precision X-Y translating stage was used to achieve finite 
fringe lines. A schematic of the experimental arrangement 
is shown in Figure 1. Because the experiment was to be a 
double-exposure hologram, the quality of the optical 
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components was not of concern. An optical bench for iso- 
lating extraneous vibrations from the experiment was also 
deemed unnecessary due to the high speed exposures being 
used. The photographic emulsions used were AGFA-GEVAERT 
10E75 holographic plates. 

2 . Michelson Interferometer 

A Michelson interferometer [Refs. 2, 6, and 8] was 
constructed to make an initial start in determining the 
feasibility and the associated problems of holography 
employing fiber optics. One path length of the interferometer 
was replaced by a fiber light guide. [Fig. 2] The refractive 
index for the optic bundle was taken as 1.5 and an appropriate 
compensation was made in the optic path lengths. The 
interferometer is a test for temporal coherence of light. 
Temporal coherence for laser sources of light is extremely 
good and is the basic quality allowing for the development 
of the modern expansion of holography. When the paths are 
exactly matched, a well-defined circular interference 
pattern should be visible on the screen. If one path length 
is changed (AL) the fringe pattern becomes less distinct. 

The point at which the fringe could be considered no longer 
existent is a subjective quality and the value of 2AL at 
this arbitrariliy accepted visibility of the fringe lines 
defines the temporal coherence length of the laser. Pulsed 
lasers have coherence lengths on the order of 10 to 20 cm. 
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3 . Holographic Arrangement 

The holographic arrangement is illustrated in 
Figure 1. The pulsed laser was mounted on an optical rail 
along one side fo the work table. In its path a removable 
front surface mirror was positioned which directed the 
alignment laser along the identical path as the pulsed laser. 
The high intensity pulse was spread by a concave lens to 
prevent the high intensity beam from scoring either the 
beam splitter or the fiber optics. The beam splitter was 
20$ reflective allowing the greater amount of energy to 
enter the fiber optics. It was necessary to direct a high 
percentage of the energy into the fibers for there is a 
^3$ internal optical loss associated with the 2k inch length. 
Also, th e, divergence of the beam leaving the fibers is 
sufficient to cause a rapid loss in beam intensity as the 
free-air path length of the beam is increased. The reflected 
beam was used as the reference beam. It was directed to a 
front surface mirror and through two concave lenses to the 
photographic surface. The other portion of the beam 
entered the fiber bundle after passing through the beam 
splitter and was then directed to the photographic plate. 

Both beams had a 30$ incidence angle to the normal of the 
photographic plate. A precision X-Y translation stage was 
used in the scene beam. It was initially used to translate 
the fiber ends and later to translate a diffusing plate 
placed in the scene beam. 
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Alignment and the relative intensities of the scene 
and reference beams were achieved with the c.w. laser. 
Alignment was achieved in the following manner. At two 
points 0.5 meters apart along the optical, rail blackened 
photographic papers were positioned. These were burned 
with the pulsed laser using a 1.5 mm. aperture. Pin holes 
were made in the center of each spot. A removable mirror 
was placed on the rail which reflected the c.w. laser beam 
through the two pin holes. This beam was a visual equivalent 
of the invisible pulsed laser beam and could be readily aimed 
through the optical components. Each beam was centered on a 
target grid set in the photographic plate holder. 

The loss of intensity through the fiber optics was 
not a controllable factor and as such, all adjustments for 
the relative intensity of the two beams were made to the 
reference beam. Two concave lenses closely spaced gave the 
needed control of the divergence of this beam. - White sheets 
of paper were held in front of each beam and a visual esti- 
mation of the intensities was made. It was attempted to 
have the reference beam from two to three times as intense 
as the scene beam. 

The equivalent free-air path lengths for both beams 
was approximately 190 cm. The free air portion of the scene 
beam, that portion from the exit of the fiber optics to the 
photographic plate, was determined by two factors; 
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having the fibers far enough away to achieve maximum 
coverage of the photographic plate, and not being so far as 
to seriously diminish the intensity of the beam. 

4 . Testing Procedure ' 

To initiate testing, all alignment devices were 
removed and a 2.5 mm* aperture was placed in the transverse 
mode selector. It proved prudent to make an additional test 
for correct alignment and relative intensities before 
proceeding with the actual experiment. This was accomplished 
by making two initial exposures, allowing only one beam at a 
time to expose each of two different plates. After develop- 
ing these plates, a more reliable visual comparison could be 
made of the two beam intensities. This also gave assurance 
that both beams were sufficiently centered on the photographic 
plate . 

To prove the feasibility of holographic interferometry 
employing fiber optics, an initial step was to. achieve 
finite fringe lines of a no-flow condition. Since the fibers 
consist of a large number of independent light sources, each 
one emanating from the end of each individual fiber, it was 
felt at first that this property of the light guide would 
obviate the need for a diffusing plate in the scene beam as 
is customarily used. As such, the fiber bundle was attached 
to the X-Y translation stage. Between exposures the stage 
was translated 0.003 inches horizontally to produce vertical 
fringe lines. Additionally, the experiment was performed 



11 



using a diffusing plate translated in the same manner, in 
place of translation of the fiber optics. 

C. RESULTS 

1. Michelson Interferometer 

The Michelson interferometric experiment produced 
positive results. A fringe pattern was observed on the 
screen from the interferometer employing the fiber light 
guide. The temporal coherence was preserved through the 
fiber optics. It became evident that the experiment was 
highly sensitive to vibration. This portion of work was 
accomplished on a vibration isolating optical bench. When 
the experiment was initially performed without fiber optics, 
the fringe pattern would waver whenever slight shocks were 
induced in the concrete floor of the .test room. Once the 
fibers were positioned, the sensitivity to vibration was 
greatly heightened. The fringe pattern was observed only 
after extensive and minute adjustments were made to the 
optical components. The fringe pattern was seen for but 
an instant. Because of the difficulty of stabilizing the 
optical arrangement, no further investigation was made to 
determine whether there was a degradation of coherency 
length by the fiber optics. The apparent increased sensi- 
tivity to vibration of the fringe formation necessitated 
using high-speed exposures achieved with a pulsed laser 
for all further work. 
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2 . Holographic Experiment 



The first phase in the experiment sought to ascer- 
tain the need for a diffusing plate. In view of the large 
number of individual fibers contained in the light guide, 
would each fiber act as an individual light source and 
produce a similar effect to that of a diffusing plate? 
Initial holograms were made of a transparent phase object, 
a rhombohedral plastic form, as well as finite fringe lines 
The technique proved unsuccessful. Finite fringe lines 
were not produced. However, a hologram of the phase object 
was achieved. It was of poor quality and contained only a 
picture of the light emitted by the fibers and several high 
intensity reflections of light from the inner walls of the 
rhombohedron. The appearance was exactly, that of previous 
holograms attempted where a diffusing plate was not 
incorporated. 

Further testing was limited to achieving a finite 
fringe pattern using the diffusing plate attached to the 
translating stage. Vertical fringe lines of a no-flow 
condition were achieved. The fringe lines were clearly 
visible, but were confined to only a portion of the 
photographic plate. 

During the alignment portion of the experiment, it 
was noted that the uniformity of the beam from the fiber 
optics was greatly influenced by the positioning of the 
incidence beam. The most uniform pattern was achieved by 
an incident beam covering approximately one third of the 
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area of the fiber ends. This result was enhanced if the 
beam was directed off center. Minute changes in its 
position on the fiber ends greatly influenced the quality 
and focus of the emitting beam. Since the pulsed beam is 
unobserved, no positive control of the beam from the fibers 
was assured. It would be reasonable to suspect this 
uncertainty for the number of failures and the eventual 
poor quality of the holograms achieved. 



II. THEORETICAL INVESTIGATION OF LIMITED FIELDS 
OF VIEW IN HOLOGRAPHIC INTERFEROMETRY 



A. INTRODUCTION 

An attempt was made to further investigate the problem 
of data reduction for supersonic flows containing opaque 
bodies. An initial effort toward this was presented by 
Jagota and Collins [Ref. 5]. Known values of density along 
various lines of sight were obtained from AGARD [Ref. 9] for 
a specific cross section in the flow about a cone. To 
diminish the effect of the presence of the cone in the flow, 
a fictitious constant density distribution equal in value 
to the density at the cone surface was added. The new den- 
sity field, which now covered the entire test section area, 
was used to calculate a fringe number array that would have 
resulted from such a density distribution. These data were 
then reinverted to obtain a density distribution which was 
compared to the original density. Reference 8 contains the 
computer program used for the data reduction. 

Further research into flow around opaque bodies has 
been facilitated by the development of an alternate computer 
program using a Fourier transform method in place of the 
Fourier expansion used in Ref. 8 to solve the pertinent 
equations. Van Houten discussed this method in Ref. 10. 

The main advantages of the program are in the computational 
time required and the numerous options available to the 
user. 
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Computational times on the order of two minutes made 
the program by Van Houten ideally suited for exploring 
different schemes to calculate the density field around 
an opaque object. 



B. PROCEDURE 

Because one would normally be starting with experimen- 
tally observed fringe data, manipulation of this information 
was chosen as a starting point instead of the density field 
change used in Ref. A fictitious density distribution 

was formulated corresponding to Pig. 8 of Ref. 10. The 
relative refractive index function was given by: 



f (X,Y) = (1 - Y 2 ) e~ 2x 



R < 1.5 



= 0 



R > 1.5 



2 2 
r = yr + y 



eq. (1) 



The computer program used this input to calculate a 
refractive index field over a ^ in. by 4 in. cross-sectional 
area. The field was calculated at 101 equispaced points on 
both the X and Y axes. Mode 1 of the program starts with 
the relative refractive index function, calculates an initial 
refractive index field at selected points from this, com- 
putes the fringe shift that would be observed from the field, 
and recalculates a new refractive index field from the fringe 
shift information. These calculations were used as standards 
against which other flow fields were to be compared. 
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Two experimental phases were studied; Van Houten's 
Fig. 8 was centered within the test field, and the same 
figure was translated 0.25 in. along both axes. The KSYM 
argument of the computer program, which is a measure of 
symmetry of the field about the origin, changed between 
cases. The translation destroyed the symmetry found along 
both axes in the original figure. 

Mode III of the program uses fringe information supplied 
as data, and from it calculates a refractive index field. 

To simulate the presence of an opaque body immersed in the 
flow, fringe information obtained from Mode I was blocked 
over a portion of the field. The blocked portion was the 
center portion of the flow field. This portion of the 
fringe data blocked by the presumed body was restored using 
two numeric schemes. The reconstructed fringe shift infor- 
mation was used to recalculate the refractive index field. 
This allowed for a direct comparison between density 
distributions of the free flow with that of a flow around 
the opaque body. 

The program from Ref. 10 was modified for this use and 
is presented in Appendix A. Two schemes were used for 
reconstructing the lost fringe information. The first 
method used a U. S. Naval Postgraduate School Computer 
Library subroutine titled INTRPL. This was incorporated in 
a subprogram [Appendix B] that read in the entire fringe 
shift information and automatically destroyed the middle 
quarter of the data. INTRPL used the remaining field* to 
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recalculate the missing portion. The second scheme [Appendix 
B] filled the same missing information with a linear inter- 
polation. The two newly constructed fringe number arrays 
were reinserted in the main computer program under Mode III 
and the resulting refractive index fields were obtained. 

C . RESULTS 

Six separate computer runs were made, the first three 
dealing with the function centered on the X- and Y-axes, 
and the last three dealing with the displaced function. 

Each group of three runs commenced with a given refractive 
index function from which was calculated the refractive 
index field at selected grid points using the Mode 1 opera- 
tion. All subsequent runs were done under Mode 3 . The 
second run of each group used the fringe shift array as 
modified by the INTRPL subprogram. The third run used the 
same information as modified by the linear interpolation 
subprogram. A diagram of the calculating procedure is 
given in Figure 3 . 

Contour plots of all refractive index fields were used 
to compare results for the different cases. Contour lines 
were produced at 0.1 intervals from the maximum function 
value of 1.0 to the minimum value. Comparisons among the 
contour plots were made using two different criteria; 
errors of the test case values from the values of the 
standard case in percent of the maximum value of the 
function, and the number of calculation grid points between 
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equal contour levels of each graph. Comparisons were 
terminated at the boundary of the simulated body. The 
body consisted of a circular region of a 0.5 in. radius 
located at the center of the cross-sectional area. 

1. Centered Function 

Figure 4 presents the refractive index field 
calculated by Mode 1 using the density function given in 
Equation 1. This is the standard of comparison for the 
other centered cases. In Fig. 5, the center portion of 
the fringe array has been removed and replaced by the sub- 
routine INTRPL. In Fig. 6, the center portion of the 
function has been removed and replaced by a linear 
interpolation. 

The greatest error encountered in comparing the 
standard, Fig. 4, with Fig. 5 was at the surface of the 
body. This error was on the order of 8$. Errors around 
most of the body were approximately 5$. Elsewhere in the 
flow, the maximum, error observed was 2.5$. Generally, the 
lower-valued contour levels away from the body were very 
closely matched. 

A comparison of Fig. 4 and Fig. 6 produced greater 
errors. The contour lines developed sinusoidal-like oscil- 
lations over the entire flow field. For comparing these 
irregular contour levels a mean value between peaks in the 
oscillation was estimated and compared with the standard 
case. Again, the greatest error was found at the body 
surface. There was a 10$ error in the positive Y-direction 
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that diminished rapidly around the rest of the body. 

Errors away from the body were negligible along the Y-axis, 
but were 6$ along the X-axis. 

2 . Displaced Function 

Figure 7, calculated in the same manner as Fig. 

3, is the standard of comparison for the displaced cases. 

In Fig. 8 and Fig. 9> the center portions of the fringe 
array were, as before.- with the centered functions, replaced 
by INTRPL and by the linear interpolation, respectively. 

A comparison of Fig. 7 with Fig. 8 showed slightly 
increasing errors over those of the centered cases. The 
greatest error occurred again where contour levels inter- 
sected the body. Errors in this region were from 7 % to 
12$, the higher error being associated with the higher 
contour levels. Elsewhere, there was good agreement between 
the two figures with the maximum error being on the order 
of 4$. 

Figure 9 again produced sinusoidal-like oscillations. 
Using the same procedure discussed for Fig. 6, the maximum 
error was found to be about 5 %* Again, as the contour lines 
neared the body, accuracy was lost. Errors had improved 
slightly in this region over the last comparison and were 
from 7$ to 10$. 

3. General Trends 

Errors for the centered functions were found to be 
concentrated along the X-axis. Errors for the displaced 
functions were generally found in quadrant III. Flow fields 



20 



for both functions showed little error along the Y-axis. It 
should be noted that the Y-axis contained a discontinuity 
and a region of negative values. 

The second method used to compare flows determined 
the number of grid points between like values of contrasted 
flows. The accuracy in the general location of contours 
from figure to figure was considered good. Contours 
locations differed by no more than two grid points. 
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III. CONCLUSIONS AND RECOMMENDATIONS 



It has been demonstrated that fiber light guides can be 
applied to holographic interferometry. The experiment was 
rudimentary and requires refinement. The problems 
encountered in the course of the experiment involved 
accurate alignment of the pulsed laser beam and attenuation 
of energy associated with the fiber optics. 

If more realistic experiments are to be conducted, 
longer fibers will possibly have to be used. Transmission 
losses of the order of 65 $ are associated with six-foot 
lengths of fibers. A means of coping with this high 
magnitude loss is needed. 

The incomplete fringe information associated with the 
limited fields of view caused by an opaque body positioned 
in a flow field was tested using two methods of calculating 
data values for the blocked portion of the fringe array. 
Both methods provided fair estimates of the flows in the 
region near the surface of the body, but there was a 
general degradation of flow information over the entire 
flow region producing errors on the order of 5 % to 10$. 

In using the Fourier transformation method for fringe 
conversion given in Ref. 10, there is a degradation to the 
flow field due to the finite spacing of the grid points. 
This effect has not fully been investigated. Errors 
associated with the two schemes for creating fringe data 



22 



for blocked portions of the flow can compound the inherent 
errors of the basic program. 

Further endeavors could be made in determining the 
effect on the investigation caused by differing the 
percentage of view that is blocked. The program could be 
run using a finer grid pattern or a greater number of views. 
These factors could all have an effect on flow errors. 
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FIGURE I. SCHEMATIC DRAWING OF THE HOLOGRAPHIC ARRANGEMENT 
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F IGURE 2. MICHELSON INTERFEROMETER WITH FIBER OPTICS 
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FIGURE 3. CALCULATION FLOW DIAGRAM 




Figure k . Standard Case for the Centered Function 
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Figure 5* Centered Function Modified by INTRPL 
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Figure 6. Centered Function Modified by INTRPL (LINEAR) 



t 



V 

v 



29 




30 




Figure 8. Displaced Function Modified by INTRPL 
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Figure 9. Displaced Function Modified by INTRPL (LINEAR) 
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